This document walks through the Reorienting to Recovery’s Recovery Definition and evaluates how the Tortoise Balanced Scenario Performs. It provides summary statistics and visualizations to evaluate recovery according to defined thresholds recovery metrics (Threshold Recovery Metrics table below) and to explore results associated with non thereshold biological metrics (Additional Biological Metrics table below).
The Reorienting to Recovery selected a few indicator populations as a starting point for evaluating recovery of scenarios. For the Fall Run Chinook model, we selected:
Upper Sacramento River
American River (?)
Stanislaus River
Tuolumne River
Indicator rivers must meet all recovery metrics in last 15 years of the simulation.
Secondarily, all independent populations must meet recovery in at least 5 years and the percentage of ind populations that meet recovery objects must be above 80% in last 15 years.
The table below summarizes Recovery Metrics that have associated Recovery Thresholds. For each of these metrics, scenarios are evaluated to see what percentage or tributaries or years they achieve the recovery target for each scenario.
| Value | Performancemeasure | Description | Target/Recovery Threshold |
| 1: Salmonid abundance | 1: Adult abundance | Represents an interest in increasing salmonid abundance across the Central Valley. Can be estimated with models as annual number of adult spawning salmonids in each tributary and ocean, tracked for 20 years. Adult abundance for each year and tributary is a direct output from the SIT DSMs and is summarized to show average Central Valley annual abundance. | 500 or more spawners |
| 2: Salmonid productivity | 2.1: Cohort Replacement | Represents an interest in having growing or stable salmonid populations. CRR can be estimated with models as the adult natural origin return ratio in each tributary, tracked for 20 years. CRR is calculated as the number of natural origin adult returns in year X+3 divided by the number of in-river spawners in year X. CRR is calculated for each year and tributary and then averaged across tributaries and years. | Mean CRR > 1 |
| 2: Salmonid productivity | 2.2: Population growth rate | Represents an alternative way to capture growing or stable salmonid populations. Can be estimated with models as average trend in population growth over a 20-year period. | Mean Growth Rate > 1 |
| 4: Salmonid genetic diversity | 4: pHOS | Represents an interest in preserving genetic diversity of natural populations and minimizing potential negative demographic effects from hatchery origin fish. Salmonid models are currently being modified to estimate proportion of hatchery origin spawners (pHOS). | pHOS < 5% |
The table below provides additional recovery metrics to help balance and refine scenarios. This set of metrics does not have clear threshold values that indicate recovery but instead provides additional context to help select more balanced scenarios.
| Value | Performance measure | Description |
| 3: Salmonid spatial structure | 3.1: Independent populations | Represents an interest in increasing and maintaining the distribution of salmonids across historically occupied watersheds. Watershed-specific abundance estimates from models can be used to calculate the number of independent viable populations in each diversity group per ESU/run. |
| 3: Salmonid spatial structure | 3.2: % of potential independent populations | Represents an alternative way to capture the distribution of salmonid, calculated as the % of potential independent viable populations in each diversity group per ESU/run. |
| 3: Salmonid spatial structure | 3.3: Dependent populations | Represents an interest in increasing and maintaining the distribution of salmonids across historically occupied watersheds. Watershed-specific CV Salmon Recovery abundance estimates from models can be used to calculate the number of dependent populations in each diversity group per ESU/run. |
| 5: Salmonid life history diversity | 5.1: Age distribution of spawning adults | Represents an interest in preserving life history diversity (age of spawners) in natural populations. Salmonid models can be used to estimate the minimum % of each age class of adults and check if targets are met (Age 4 >35%, Age 5+ >20%). |
| 5: Salmonid life history diversity | 5.2: Size distribution of juveniles | Represents an interest in preserving life history diversity (size distribution of juveniles) in natural populations. Salmonid models can be used to estimate the variation in juvenile abundance of each life stage (fry, parr, yearling) and calculate variation across years. Size distribution of juveniles is summarized in a single metric by using a Shannon diversity index. |
| 5: Salmonid life history diversity | 5.3: Size and timing of juveniles at ocean entry | Represents an interest in preserving life history diversity (size and timing of juveniles at ocean entry) in natural populations. Size distribution and timing of juveniles at ocean entry is summarized in a single metric by using a Shannon diversity index. |
| 5: Salmonid life history diversity | 5.4: Habitat diversity | Represents an interest in preserving habitat diversity that facilitates life history diversity in natural populations. Could be estimated as the amount and relative % of available habitat of different types (measured in area and days). Specific methods using the salmon models to calculate and compare optimal habitat ratios to population growth are being developed. |
| 6: Abundance, productivity, and distribution at levels that support ecosystem health | 6: Ecosystem health | Represents an interest in maintaining resilient and functioning ecosystems due to adequate salmonid abundance. Could be captured with a proxy metric of ratio of spawner abundance to habitat (marine derived nutrient proxy). The marine derived nutrient proxy metric gives the grams of marine nutrient per square meters of river. To calculate marine derived nutrient proxy, we calculated the total acres by multiplying habitat extent lengths by river width. We then multiplied the average annual number of spawners by the average weight (assumed 21 grams) and divided by our total area CV Salmon Recovery Draft for Planning Team Review 11 (square meters) to get a grams of marine derived nutrient per square meters. |
| 7: Biological recovery objectives met by certain time | 7: Time to biological recovery | Represents an interest to recover salmonids quickly. Can be estimated by salmon models as the # of years until recovery objectives are met. See Phase 1 report for a list of recovery objectives. |
The Tortoise scenario includes efforts in each of the 4 knobs that we can turn in the R2R process: Habitat, Hydrology, Harvest (oceans & rivers) and Hatcheries. The table below summarizes the Tortoise scenario in comparison to the Baseline scenario.
| Recovery scenario | Habitat | Hydrology | Harvest (oceans & rivers) | Hatcheries |
| Baseline | Current + planned | Current operations | Current levels | Current operations |
| Tortoise | Current + planned Floodplains / rice fields (Sac/SJ) Food subsidies (dry years) Predation reduction (large-scale, dry years) |
Ecological Functional Flows (EFF) Dry years |
No harvest of dry year cohorts (ocean, in-river) Harvest of hatchery fish only (ocean, in-river) Tribal harvest prioritized All years |
Phased hatchery Weirs All years |
The following code block runs the Tortoise scenario.
# Tortoise --
new_params_tortoise <- R2Rscenario::load_scenario(R2Rscenario::scenarios$balanced_scenarios$tortoise,
species = "fr")
new_params_tortoise$movement_hypo_weights <- c(1, rep(0, 7))
# BASELINE --
new_params <- fallRunDSM::r_to_r_baseline_params
new_params$movement_hypo_weights <- c(1, rep(0, 7))
new_params$san_joaquin_flows <- matrix(0, nrow = 12, ncol = 21,
dimnames = list(month.abb, 1980:2000))
# seed
r2r_seeds_tortoise <- fallRunDSM::fall_run_model(mode = "seed",
seeds = fallRunDSM::adult_seeds,
..params = new_params, # seed with baseline, then Tortoise
delta_surv_inflation = FALSE)
# run model
r2r_model_results_tortoise <- fallRunDSM::fall_run_model(mode = "simulate",
scenario = "tortoise",
..params = new_params_tortoise,
seeds = r2r_seeds_tortoise,
delta_surv_inflation = TRUE)
Model results should be evaluated relative to the Baseline model. The following plots show the % difference from baseline results for the Tortoise scenario. You can look at the tab to the right to see the total spawner trends of the Tortoise Scenario in comparison to the Baseline Scenario. (Make add a percent difference table in here too) - Think through the best way to present these. Also add a percent difference bar plot)
The two tabs below show plots that describe if a tributary meets the recovery objective.
Abundance and Genetic Diversity metrics are evaluated for each Independent population in every year of the scenario.
Productivity Metrics are averaged for each independent population over the last 10 years of the simulation. Productivity metrics must on average for each independent population, meet recovery thresholds to consider the system at recovery.
(TODO: maybe color indicator pops differently)
## Rows: 1,240
## Columns: 12
## $ location <chr> "American River", "American River", "America…
## $ year <dbl> 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 1…
## $ scenario <chr> "Tortoise", "Tortoise", "Tortoise", "Tortois…
## $ run <chr> "fall", "fall", "fall", "fall", "fall", "fal…
## $ performance_metric <chr> "1 All Spawners", "1 All Spawners", "1 All S…
## $ value <dbl> 56460, 37538, 124134, 244500, 169989, 511939…
## $ origin <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ size_or_age <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ month <int> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
## $ `Spawners Present` <dbl> 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,…
## $ diversity_group <fct> 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,…
## $ div_group_with_watersheds <chr> "5:", "5:", "5:", "5:", "5:", "5:", "5:", "5…
## # A tibble: 4 × 1
## location
## <chr>
## 1 Battle Creek
## 2 Bear Creek
## 3 Cow Creek
## 4 Upper Sacramento River